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Abstract 

We study an assembly line balancing problem that occurs in shel- 
tered worker centers for the disabled, where workers with very different 
characteristics are present. We are interested in the situation in which 
parallel assembly lines are allowed and name the resulting problem as 
parallel assembly line worker assignment and balancing problem. We 
present a linear mixed-integer formulation and a four-stage heuristic 
algorithm. Computational results with a large set of instances recently 
proposed in the literature show the advantages of allowing alternative 
line layouts. 

Keywords: parallel assembly line balancing; heterogeneous workers; heuris- 
tics. 



1 Introduction 

The current economic downturn raises substantial challenges to the goals 
and actions proposed for the Disabled integration. According to the Inter- 
national Labour Organization (ILO), people with disabilities represent an 
estimated 10 per cent of the world's population, where almost 500 million 
are of working age. Labour market activation for their inclusion is easier to 
address in periods of increasing labour demand than in times of recession. 
Nevertheless, there is evidence of national actions to create more flexible 
work solutions, and the success of such initiatives show that there is scope to 
learn from best-practice cases, including a comprehensible knowledge about 



how and why they work (Greve, 2009) 



In this sense, one of the actions most commonly adopted to facilitate the 
integration of people with disabilities into the labor market has been the 
creation of Sheltered Work centers for Disabled (henceforth SWDs). This 
model of socio-labor integration tries to move away from the traditional 
stereotype that considers these people unable to develop continuous profes- 



sional work (Miralles et al. 2010). Just as any other firm, a SWD competes 
in real markets and must be flexible and efficient enough to adapt to market 
fluctuations; the only differences being that the SWD receive some govern- 
mental help to compensate the fact that most of its workers (usually around 
70%) are disabled, and some limitations and therapy requirements must be 
taken into account. 

In countries such as Spain, this labor integration formula keeps being 
successful in offering jobs to disabled people. One of the strategies of SWDs 



to facilitate this integration is to use assembly lines as the most accessible 



configuration. In this sense Miralles et al. (2007) were the first to evidence 



how the division of work in single tasks enables many possible job assign- 
ments that make invisible the disabilities. 



1.1 Literature review 

Traditional assembly line balancing research has focused on the simple as- 
sembly line balancing problem (SALBP), that uses several well-known sim- 
plifying hypotheses, which reduce the complex problem of assembly line 
configuration to the "core" problem of assigning tasks to stations so that 
certain precedence constraints are fulfilled. 

In recent years a growing amount of literature has been produced dealing 
with more realistic and generalized problems. Regardless of the theoretical 
importance of the SALBP, balancing real- world assembly lines requires the 
observation of a large variety of additional technical and organizational as- 
pects, which can heavily affect the structure of the problem. In fact, in the 
last decade a big effort has been made towards modeling real world assem- 
bly line systems through different extensions of SALBP, aiming to narrow 



the former gap between research and practice (see reviews of Scholl and 



recently Dolgui and Battaia (2013)). 



Becker (2006) , Becker and Scholl (2006), Boysen et al. (2007, 2008) or more 



In this sense the so-called assembly line worker assignment and balanc- 
ing problem (ALWABP) represents one of these recent efforts made by the 
Academia. The ALWABP focus on the heterogeneity of task times and 
the presence of incompatibilities, defining a new set of realistic hypotheses 
inspired by the SWD assembly lines where disabled workers usually exe- 
cute tasks at diflFerent rates (iMiralles et al.l 120071). The ALWABP defines 



worker-dependent processing times, which allow taking into account the lim- 
itations of each worker and can, therefore, be useful in environments other 
than SWDs where workers also have diverse speeds to perform certain tasks 



(Araiijoet al. 2012) 



Since the initial paper of Miralles et al. (2007), many other references 



have contributed to give this problem visibility in the scientific literature and 
several solution methods have been proposed. Exact methods have focused 



on branch-and-bound strategies (Miralles et al. 2008; Borba and Ritt 2012) 



although most of the efforts have concentrated on heuristic and metaheuris- 
tic methods due to the NP-hard condition of the problem, using strategies 



such as clustering search (Chaves et al. ; 2009), tabu search (Moreira and 



Costa; 2009), Iterated Beam Search (Blum and Miralles 2011), construe- 



tive heuristics (Moreira et al. , 2012) or genetic algorithms (Moreira et al. 



2012; Mutluet ah 2013). 



1.2 Contribution and outline of this work 

Despite this intense research on the ALWABP, the hterature stiU presents a 
significant gap in terms of proposing flexible configurations in the context of 
SWDs (or in contexts with heterogeneous workers in general). In the par- 
ticular case of SWDs that motivates this paper, managers require easy tools 
to handle the workers heterogeneity while optimizing productive efficiency. 
Therefore, it is very important for them to count on as many alternative 
combinations as possible for job assignments with the workers available. 

In spite of these flexibility requirements, the ALWABP is traditionally 
limited by the hypothesis of a single serial assembly line that maximizes 
the production rate with the resources (workers) available. We propose to 
open this scope and analyze alternative configurations that might adequate 
better to the available workforce. In particular, we study the possibility 
of assigning workers to teams, which can then operate parallel assembly 
lines. Our hypothesis is that by exploring these new configurations, a large 
number of new worker-tasks assignments will appear and possibly allow the 
obtention of higher productivity levels. 



To the best of our knowledge, Araujo et al. (2012) were the sole authors 



to analyze such possibility. The authors study the assignment of workers in 
parallel, at the same station of a single line. Although some positive results 
have been obtained, this approach has the drawback of possibly increasing 
the complexity of coordination and control. In this paper, we propose to cre- 
ate new (parallel) assembly lines which will make use of all available workers. 
This is done by building work teams which can operate independently from 
each other. 

In this sense, this paper will completely define an extension of ALWABP 
that can provide better solutions than the traditional approach, just by 
allowing parallel complete lines. We name this extension as the parallel 
assembly line worker assignment and balancing problem (PALWABP). In 
the following, this problem is analyzed by means of an example. Then, a 
formal definition and a mathematical model are proposed. In Section |4j a 
heuristic methodology for obtaining good-quality results in low computation 
times is detailed. Computational results over a large set of instances are 
presented and analyzed in Section [5| Finally, general conclusions and hints 
for future works end this paper in Section [6} 



2 Parallel Assembly line Worker Assignment and 
Balancing Problem 

As introduced in section [T| the ALWABP is a generalization of the SALBP 
where, in addition to the assignment of tasks to stations, a set of hetero- 
geneous workers that execute the tasks also has to be assigned to stations 



dMiralles et al.t|2008D . 

Most previous proposals in the literature face the ALWABP of type 2 
(minimization of cycle time given a fixed number of workers) , the most typ- 
ical situation in reality, and have been mostly evaluated with the set of 



320 benchmark instances first proposed by Chaves et al. (2009). In order 
to illustrate the extension here proposed, we have extracted the problem 
HESKIA_64 of this benchmark and expressed its input data and the opti- 
mal solution in Table 1 and Figure 1, respectively. In the table, for every 
task (rows) several operation times are possible depending on the worker 
(columns). If a task is considered unfeasible for certain worker, the incom- 
patibility is represented by a dash in the corresponding cell of the input 
data matrix. In the optimal solution of Figure [T] we can notice how all the 
precedence constraints are respected while ensuring an optimal cycle time 
of 126 s (given by the bottleneck station, which is the second station with 
worker W6 performing task 20): 

This is the minimum cycle time we can get with the traditional ALWABP 
approach. But if we allow the possibility of designing complete parallel lines 
with the available workers, then solution space is enlarged, and lower global 
cycle times can be reached. To illustrate this fact, let's focus now on table 
[I] (but now taking into account the light shaded/dark shaded cells, and the 
bottom of the table): 

In the table we have represented an alternative solution by using a light 
shade on cells with task- assignments within line 1, and a darker shade on 
those cells with task-assignments within line 2; and where it can be easily 
checked how precedence constraints have been respected in both parallel 
lines. At the bottom of the table we summarized the workstations of each 
sub-line with the following independent work teams: 

• Line 1 is composed of four consecutive stations with workers W5, W4, 
Wl and W7. Workers Wl and W4 are bottlenecks with a cycle time 
of 135. 

• Line 2 is composed of three consecutive stations with workers W3, W2, 
and W6. Workers W2 and W6 are bottlenecks with a cycle time of 
354. 




Figure 1: Optimal solution for HESKIA_64 with the traditional ALWABP 
approach 



To compute the global cycle time (CT) of this configuration we can define 
the corresponding Throughput Rates (TR) of both lines. In this case, we 
have: 

• Line 1 has CTi = 135 s/part, which is equivalent to a throughput TRi 
= 26.67 parts/h 

• Line 2 has CT2 = 354 s/part, which is equivalent to a throughput TR2 
= 10.17 parts/h 

Thus, the global Throughput Rate using this parallelized configuration 
is: TR = TRi +TR2 = 36.84 parts/h, that corresponds to a combined cycle 
time CT = 97.7 s/part 

Therefore, the optimal cycle time of 126 obtained for HESKIA_64 with 
the traditional serial ALWABP approach is improved by more than 22% 
just by assigning workers to parallel complete lines. As will be later demon- 
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Table 1: Matrix of process times for HESKIA_64 and solution with two 
parallel lines 



strated, this is not an isolated successful example, since many others bench- 
mark cases get better cycle times through this strategy. 

In addition to this improvement in global efficiency, the proposed par- 
allel assembly lines layout has also some practical advantages, that can be 
summarized as follows: 

• In the case of workers with very slow operation times, it is much easier 
to integrate them if we allow parallel lines. 



With the parallehzation, the facihties and tools must be replicated, 
with the possible inconvenient of certain physical constraints or instal- 
lation costs but with a clear advantage: the manager gets many more 
possible task assignments, that can be crucial for job rotation purposes 



(see Costa and Miralles (2009) and Moreira and Costa (2012)) 



In fact, having more lines supposes that workers execute a higher num- 
ber of tasks, in average, avoiding work routine and improving workers 
capabilities. 

As it happens with SALBP, when we allow parallehzation the cycle 
time can be lower than the largest operation time. In the case of 
ALWABP, lower than the slowest task to execute (obtained with a 
max-min operation over all worker x tasks execution times). 

Furthermore, parallelizing complete lines, instead of single stations 
and/or tasks, makes the assembly line management and self-control 
easier: every operator is assigned to a work team in charge of certain 
station within a complete parallel serial line, with independent feeders, 
work in progress stocks, and tools. 



2.1 Review on parallelization within assembly lines 

Despite the fact that most assembly lines addressed in the literature are quite 
different from ALWABP, it is important to review the most important refer- 
ences that somehow face parallelization. Thus, we can start chronologically 
by citing the work of Askin and Zhou (1997), who consider mixed-model 



assembly lines with task dependent equipment costs and arbitrary number 
of parallel workstations. They propose a heuristic procedure that uses a 
threshold value for the equipment utilization. Siier ( 1998 ) proposes a sim- 



ple resolution method in three phases, while McMullen and Prazier (1997 



1998 ) propose a heuristic and a simulated annealing technique for a problem 



with parallel stations and stochastic times. In both cases, they measure the 
line performance by the total cost and the extent to which cycle time is 



met. Vilarinho and Simaria (2002) introduce a two-stage heuristic method 



for mixed-model assembly lines, where their primary goal is to minimize 
the number of workstations for a given cycle time, and the secondary goal 
is to balance the workload between workstations. IBukchin and Rubinovitzl 



(2003) study the equipment selection problem on parallel workstations, in- 



vestigating the influence of assembly sequence flexibility and cycle time on 



the balancing improvement due to the station parallehng. Simaria and Vi 



larinho (2004) propose a model and a genetic algorithm for a mixed model 



assembly line with parallel stations. Two years later the same authors pro- 



pose in (Vilarinho and Simaria 2006) an ant colony optimization algorithm 



for balancing the mixed-model assembly lines with zoning restrictions and 



parallel workstations. Boysen and Fliedner ( 2008 ) present a versatile graphic 



algorithm that can be adapted for parallel stations, and Ege et al. (2009) 



consider parallel workstations and propose a branch and bound algorithm 
to minimize total equipment and workstation opening costs. 



It is important to highlight (Scholl and Boysen, 2009), that consider 



a problem first introduced by [Gokgen et al" (2006) and named as parallel 
assembly line balancing problem (PALBP) where different products and/or 
models can be assigned to different serial lines arranged in parallel while 
considering a joint cycle time. From 2009, the PALBP has been the main 
topic of research in parallelization where, even not being the focus of this 
paper (since in our case we have the same product to be performed in par- 



allel assembly lines), we can highlight two recent references: Kara et al. 



(2010) proposed a fuzzy goal programming model that can be used for bal- 



ancing this kind of multi-product parallel lines. Ozbakir et al. (2011) pro 



pose a multiple-colony ant algorithm for balancing multi-product with two 
objectives (maximizing the line efficiency and minimizing the idle time of 
workstations) parallel assembly lines; obtaining good results with respect to 
other existing algorithms. 



Finally, Kellegoz and Toklu (2012) propose efficient branch and bound 



algorithms to solve assembly lines with parallel multi-manned workstations, 
a configuration that is necessary for big-sized complex products containing 
a large number of assembly tasks, most of which have long task times. 

The only reference that studies the ALWABP with parallelization is 
(Araiijo et al. : 2012), where two extensions are explored: in the first ex- 



tension they allow multiple stations to execute the same set of tasks in the 
same stage; and in the second extension different workers can complement 
their capabilities collaborating in the same product within certain worksta- 
tion. In fact, the first extension uses parallelism to get more combinations of 
worker-tasks assignments. However, the difference with respect to our pro- 
posal is that we force to configure complete parallel assembly lines, building 
up independent work teams that don't increase complexity of control as 



it may happen in (Araiijo et al. 2012); where parallelizing only certain 



stations and/or tasks can be cluttering because the production flow is di- 
vided/duplicated (what can be hard to coordinate, especially in the case of 
mental disabilities). 



These practical insights made us reflect on the need of proposing models 
and algorithms for complete parallel stations in ALWABP. In the following 
section, we present a formal definition of this problem by means of a classical 
literature taxonomy and also via a mathematical model. 

3 Formal definition and mathematical model 



We formally classify the PALWABP using the nomenclature of Boysen et al. 



(2007), that structures the vast field of assembly line balancing problems by 



means of a notation consisting of three elements [a|/3|7], where: 

• a concerns the precedence graph characteristics; 

• /3 concerns the station and line characteristics; 

• and 7 concerns the optimization objectives. 

Various possible values are defined for each of these three elements, cov- 
ering all existing assembly line balancing problems in the literature, and 
where the absence of values indicates the assumption of SALBP hypothesis 
(e. g. SALBP-2 is coded as [||c]). 

[Boysen et al. (2007) classified the ALWABP-2 as [pa,link,cum|equip|c]. 



In the PALWABP, we allow multiple assembly lines (/3=pline). Therefore, 
the PALWABP-2 can be classified as [pa, link, cum|pline, equip |c]; whereas the 



previous proposal of Araiijo et al. (2012) has formal and practical differences 



that have been properly detailed, being coded as [pa,link,cum| pstat,equip|c]. 

3.1 Mathematical model for PALWABP 

Let (A^, <) be a partially ordered set of tasks in which i < {>)j indicates that 
a task i must precede (succeed) a task j . Also let W he a, set of workers and 
Pwi an integer associated with each pair (w, i) G (VF, A^) indicating the time 
worker w spends to execute task i. Also, /«, is the set of tasks that worker w 
is not able to perform. Consider an assignment at : N ^ S of the tasks to 
a linear sequence of stations S = {1,2,..., m} respecting the tasks partial 
order and an assignment a^, : VF — )• 5" of the workers to the same linear 
sequence of stations S = {1,2, .. . ,m}. The load of a station is the time 
the worker assigned to that station needs to execute those tasks assigned to 
the station. Finally, the cycle time of such assignments is the largest load 
among all stations. The ALWABP-2 aims at finding assignments at and a^ 
minimizing the line cycle time. Let ALWABP-2(A^, W,pwi) be such optimal 
cycle time. 

10 



The PALWABP-2 considers the same input data and an extra integer, 
Kmaxi which hmits the maximum number of parallel assembly lines. It 
then aims at finding a partition of W in Wi^ W2---Wk^^x sets, such that the 

overall throughput rate J2k=i'' alwabp 2(n w 1 ^® maximal. We assume 

that a set Wk can be empty and, in this case, ALWABP-2(iV, Wk,Pwi) = oo- 

In the following, we present a mathematical model for the PALWABP-2. 
We define the following set of variables: 

Xswik Binary variable. Equals one if worker w executes task i in 

station s of line k and zero otherwise 
yswk Binary variable. Equals one if worker w is designated to 

station s in line k and zero otherwise 
Zk Binary variable. Equals one if line k is active (that is, if at 

least one worker is assigned to that line) and zero otherwise 
Cfc The cycle time of line k 

We can write the following model for the ALWABP with parallel lines: 



Min — f^ 

E"-Max _J_ 
P=l C'k 

Subject to: 
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The objective function M minimizes the combined cycle time. Con- 
straints Q guarantee that in every active hne, each task is executed by 
a single worker in a single station. Constraints ([3| guarantee that every 
worker will be assigned to only one station from one of the parallel lines. 
Constraints Q allow a single worker in each station of active lines and zero 
in each station of non-active lines. Constraints ([5]) establish the precedence 
relations between the tasks in each assembly line. Constraints (pi) state that 
a worker can execute tasks in one stage only if that worker is assigned to 
that stage. Constraints M define the cycle time for each assembly line as 
the execution time among all workers from that line while constraints (fsl) 
handle the task- worker infeasibilities. 

The objective function (IT]) makes this model non-linear. In order to 



linearize it, we adapted the linearization proposed by Araiijo et al. (2012) for 
the ALWABP with parallel workstations. In that work, the authors changed 
the objective function to maximize the throughput rate of the assembly line, 

12 



which is equivalent to minimizing the cycle time. This can be adapted for 
the PALWABP by changing the objective function to maximize the sum of 
the throughput rates of each assembly line. Using the new variables: 

-Ffc The production rate of line k. 

fswk The production rate of worker w in stage s of line k. Equals 
zero if the worker is not assigned to that stage of that line. 

Vswik Auxiliary variable used for linearization. Corresponds to 

Jswk * -^xwik 

We can write a linear model for the PALWABP as: 
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This model is used to solve small instances (see Section pi). Neverthe- 
less, due to its complexity, a heuristic solution method is developed in the 
following section. 

4 A four-stage heuristic solution method 

In this section we present a 4-stage heuristic to solve the PALWABP-2. 
The core idea of the method is the generation of sets of workers that can be 
assigned together in an independent line. Each of such team of workers must 
be able to execute all tasks while respecting the precedence constraints. 

The heuristic procedure is divided into four main steps, as shown in 
Figure [2] First the sets of tasks each worker can execute are obtained. 
Secondly, with these sets, feasible worker sets (i.e., workers that can be 
assigned together to a complete line) are generated. In the third step, each 
set is evaluated according to the production rate of the obtained assembly 
line. Finally, a maximum of Kmax sets are selected to compose the final 
solution. In the following, these four steps are detailed. 



Feasible worker-tasks sets generation 




' 


' 






Feasible worker sets generation 








' 


' 








Sets evaluation 






' 


' 






Sets selection 





Figure 2: Flowchart for the heuristic for the PALWABP 



4.1 Feasible worker-tasks sets generation 

In this first phase, the sets of all tasks that can be executed in a single 
station by each worker are obtained. If cycle time constraints are ignored, a 
worker w can execute any set of tasks as long as he is not obliged to execute 
a task in /^„. Using binary variables Xi equals one if a task i belong to a set. 
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the following constraints express this condition: 

x^ = 0, Vi E /^, (22) 

Xi + Xj < 1, \/i,j G N\Iw, k e Iw\i < k,j > k. (23) 

Constraints ( |22[ ) explicitly forbid the selection of tasks belonging to I^ while 
constraints ( |23| ) indicate that a set cannot contain two tasks if this will 
require the execution of a task in 7^^,. 

In order to avoid sets of tasks that are contained in other sets, we can 
solve the following linear problem: 



Max^: 



(24) 



subject to (22) and (23) 



This will generate the set with maximum number of tasks. To generate 
sets with tasks not included in this solution, it suffices to re-solve the model 
with the constraint Xj = 1, at most once, for all task i G N\Iu]. 

Alternatively, we propose a simple constructive procedure. Let i £ N he 
a task and / C A^ be a set of tasks. We use the notation i '^ {iC}I to indicate 
that task i does not succeed (precede) any of the tasks of /. The method 
presented in Algorithm [T] is able to quickly determine T^„, the sets of tasks 
that can be executed by a worker w. 

Algorithm 1 Task sets generation 
Require: N,w ^W 



for all I <Z Iw do 

s = {ieN\I^\iitI,i:;^ I^\I} 
if s (^ j £ Tw then 

end if 
end for 

Return T,,, 



The algorithm receives as input the worker being analyzed and the or- 
dered set of tasks. It then considers all subsets of tasks that worker w can 
not execute (line 2). In its next step, it assumes that these tasks have been 
executed in stations previous to the station to which w has been assigned 
and adds to the set of tasks that w can execute all tasks that do not succeed 
any other task that w can not execute (line 3). All the sets of feasible sets 
of tasks are collected in set T^ (line 5) which is then returned (line 8). 
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Example: Consider the precedence graph shown in Figure |3j Worker wi 
can not execute task 2 while worker W2 can not execute task 3. 





45 

(IMZ)GMIMZ) 



4s 



worker 1 



worker 2 



worker 3 



Figure 3: Example of precedence graph 

The subsets of Iwi and I^j (see line 2 of Algorithm [l]) are {0, {2}} and 
{0,{3}}, respectively. They yield T^i = {{1, 3, 5}, {3,4, 5}} and T^^ = 
{{1, 2,4}, {2,4, 5}}. Worker w^ can execute all tasks and, therefore, T^g = 
{{1,2,3,4,5}} 



4.2 Feasible worker sets generation 

The goal of this phase is to select one set of tasks from each worker so that 
each task is covered by at least one set. By doing this, we want to ensure 
the existence of a viable solution for the ALWABP with these workers. The 
information obtained in the previous stage is now used to determine all 
combinations of workers that can be assigned together to a team. For each 
worker, let T^ be the sets of tasks obtained and aiwt be a parameter equals 
to one if task i belongs to set t £ Tw Using variables y^t (equals one if 
worker w is selected and executes subset t*'^ of his set of tasks), a feasible 
worker set must respect the following constraints. 



i=l 
w&W t=l 



yw£W, 



yi£N. 



(25) 



(26) 



Constraints (25) indicate that a worker w, if selected, can only execute tasks 



that are a subset of t G T^. Constraints (26), in turn, guarantee that all 



tasks can be executed by at least one worker in the line. 

In this situation, we are interested in obtaining the sets of workers that 



respect constraints (25) and (26) with minimum cardinality, because a line 
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remains feasible if a new worker is added to it. A minimum cardinality set 
can be obtained by solving the problem: 

IT™ I 
MinJ^J^y^i (27) 

w&W t=l 



subject to (25) and (26). 



A set of workers Ws configuring a feasible line can be obtained by simply 
choosing all workers for which Xlt"' V'^t ~ ^- ^^ before, other sets can be 
obtained by repeatedly solving this problem while excluding the already 
obtained solutions. For example, if a set with workers wi and W2 have been 
obtained, one can search for a new set solving the problem above with the 

addition of a constraint X^f=™^ Vwit + Z]t=T y«'2t — 1- 

Algorithm 2 Subproblems generation 
Require: T^^, 

1: Wr = 0, setJist = 0, worker _list = 

2: i = 

3: while 3w GW — worker list, s G r^|i G s do 

4: set-list = set-list U {s} 

5: worker-list = worker list U {w} 

6: if {{N - set-list) = 0) then 

7: Wr = WrD {workerJist} 

8: Remove last element from workerJist and set-list 

9: else 

10: Select i £ (N - set-list) 

11: end if 
12: end while 
13: if setJist / then 

14: Remove last element from workerJist and setJist 
15: goto [3] 
16: else 

17: Return Wr 
18: end if 

Alternatively, we propose Algorithm [2] to generate such sets. This algo- 
rithm enumerates all possible combinations of workers that may be assigned 
to a team. It starts with an empty list of sets of tasks. In each iteration, 
if there is a task that has not yet been executed, we select a set of tasks 
containing this task (linep|. Both the set and the worker that executes the 
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set are added to their respective lists (lines H^ and p| . If all tasks are exe- 
cuted, the subproblem is complete and is added to the list of subproblems 
to be solved in the next phase (line w\ . Observe that if there are workers 
remaining, any worker we add to the list will also result in a valid solution. 
If there are no remaining workers or if the remaining workers cannot execute 
the task, the method tries another subset of tasks of the last worker in the 
workers list. If all subsets of a worker were already tested, we remove this 
worker from the worker list and the method continues with another worker. 
This process continues until we generate all the possible combinations of 



workers (line 17). 

Observe that Algorithm [2] generates sets of workers with minimum car- 
dinality. In order to generate all feasible sets we can either suppress line [8] or 
add the remaining workers to each set in Wr- The last option is preferable 
to avoid redundant sets. 



Example (continued): For the example of Figurepl all possible obtained 
combinations of sets are: 



Line 


worker wi 


worker W2 


worker w^ 


1 


r^,(i) = {1,3,5} 


T^,(l) = {1,2,4} 




2 


r^,(i) = {1,3,5} 


T^,(2) = {2,4,5} 




3 


r^,(2) = {3,4,5} 


T^,(l) = {1,2,4} 




4 


T^,(l) = {1,3,5} 




r^3(l) = {1,2,3,4,5} 


5 


r^,(2) = {3,4,5} 




r^3(l) = {1,2,3,4,5} 


6 




T^,(l) = {1,2,4} 


r^3(l) = {1,2,3,4,5} 


7 




T^,(2) = {2,4,5} 


r^3(l) = {1,2,3,4,5} 


8 


r»,(l) = {1,3,5} 


r^,(l) = {1,2,4} 


r^3(l) = {1,2,3,4,5} 


9 


r^,(i) = {1,3,5} 


T^,(2) = {2,4,5} 


r^3(l) = {1,2,3,4,5} 


10 


T^,(2) = {3,4,5} 


r^,(l) = {1,2,4} 


r^3(l) = {1,2,3,4,5} 


11 


r^,(2) = {3,4,5} 


T^,(l) = {2,4,5} 


r^3(l) = {1,2,3,4,5} 


12 






r^3(l) = {1,2,3,4,5} 



Each of these combinations result in at least one feasible ALWABP so- 
lution. 

4.3 Sets evaluation 

The two first phases of the method generate a series of subproblems. The 
next phase of the method will evaluate these subproblems in order to find 
partial solutions. Each subproblem consists of a serial ALWABP with work- 
ers Wo C W and the solution of a subproblem is a partial solution for the 
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PALWABP-2. In this phase, these subproblems are solved in order to obtain 
the productivity of each possible line. This can be done using the model 



proposed by Miralles et al. (2007) on each subset, which uses binary vari- 
ables Xswi equals one if task i is assigned to worker w at workstation s and 
binary variables ysw equals one if worker w is assigned to workstation s. 



Alternatively, we use the fast constructive heuristic developed by Moreira 



et al. (2012). This heuristic calculates a lower bound for the cycle time and 



tries to generate a solution with this cycle time. If no feasible solution is 
found, the cycle time is increased by one unit and the procedure is repeated 
until either a feasible solution is found or an upper bound for the cycle time 
is reached. Task and worker assignments are decided sequentially based on 
heuristic priorities. 



The cycle time search used by Moreira et al. (2012) is a simple sequential 



search starting with a known lower bound and increasing one unity in the 
cycle time at each iteration until a feasible solution is found. In order to 
further improve the speed of the method, at each iteration, we increase the 
cycle time by the smallest task time among all tasks and workers, that is, 
c = c + min^ tm . This greatly improves the method speed with little effect 
on solution quality. 

Example (continued): Optimal assignments and the associated cycle 
times for each of the 12 lines obtained after step 2 are listed in the table 
below. When using the constructive heuristic, suboptimal results may be 
found. 



Line 


CT 


worker wi 


worker W2 


worker ws 


1 


10 


3,5 


1,2,4 




2 


9 


1,3 


2,4,5 




3 


10 


3,5 


1,2,4 




4 


9 


1,3 




2,4,5 


5 


10 


3,5 




1,2,4 


6 


9 




1,2,4 


3,5 


7 


10 




2,4,5 


1,3 


8 


8 


3 


1,2,4 


5 


9 


7 


1 


2,4,5 


3 


10 


7 


5 


1,2,4 


3 


11 


8 


3 


4,5 


1,2 


12 


18 






1,2,3,4,5 
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4.4 Sets selection 

After the third step, we have a set of partial solutions associated with subsets 
of workers. The last phase of the method combines these partial solutions 
to create a solution for the PALWABP-2. Let Wr = {Wi, VF2, • • • W\Wr\} 
be all sets of feasible sets of workers. Also, let TR^ be the throughput rate 
associated with a line with workers Wr and q^r be a parameter equals one 
if worker w belongs to set Wr ■ With the use of variables Cr equals one if set 
Wr is selected, the problem of partitioning the workers in a most productive 
lines scheme can be written as: 

Max ^ TRrCr (28) 

r&Wr 



subject to 



\Wr\ 

^ qyjrer = 1, Vu; G W. (29) 



r=l 



In which constraints ( 29 ) simply indicates that each worker must belong to 
a single selected set. 

Observe that, for any subset of workers Wi associated with a partial 
solution, Wi U {w} is also associated with a partial solution. Therefore, if 
{Wi, W2, ..., Wn] is a feasible solution, then {Wi, {W2 U ... U Wn}} is also a 
feasible solution. Algorithm [3] uses this idea to enumerate the combinations 
of sets of workers. We use the notation UWj to indicate the union of all sets 
in Wi. 

Algorithm 3 Sets selection 
Require: Wr 



solutions = 
for all Wi C Wr do 
if UWi = W then 

if Wi n W2 = VWi, W2 G Wi then 

solutions = solutions U Wi 
end if 
end if 
end for 
Return solutions 



The algorithm analyzes all combinations of sets of workers and checks 
if each worker is present in exactly one assembly line. Each combination of 
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sets that satisfies this condition is a solution for the PALWABP. The final 
phase of the method will select the solution with the lowest cycle time. 

Example: The best serial lines obtained (using the three workers) are lines 
9 and 10, with a cycle time of 7. Nevertheless, by using worker wi and W2 in 
line 2 and using worker w^ in line 12 (alone), a cycle time of CT = -j — j- = 6 

9 + T8 

is obtained in this toy example. 

5 Computational results 

In this section we present the results obtained using the mathematical model 
and the heuristic method described in sections [3] and |4j respectively. To test 
these strategies, we generated a set of instances for the ALWABP based on 



the instances provided by Otto et al. (2013) for the SALBP. These instances 



vary in graph structure (bottleneck, chain or mixed), 'trickiness', order of 
strength and times distribution; representing a much more robust SALBP 



benchmark than the classical benchmark of Hoffmann (1990), as demon- 
strated by the authors. For each instance, we generated four instances for 
the ALWABP, with different time factors (2 or 5) and infeasibility rates (10% 
or 20%). A time factor of n means that, for every task with execution time 
of t in the original instance, the execution time of each worker is a random 
integer between t and n * t. The number of workers in each instance is equal 
to the upper bound on the number of stations in the original instance. We 
used the sets of instances with 20 tasks and 0.2 order strength as our base 
for generating a total of 900 instances. With this we have a more complete 



benchmark than the one provided by Chaves et al. (2009), allowing a better 



analysis of methods and extensions of the ALWABP. 

Both strategies were implemented and tested in a computer with a Intel 
Core 2 Duo T5450 processor, 1,66 GHz and 3 GB of RAM. The following 
sections detail the results obtained. These are compared with the best know 
solutions for the respective serial ALWABP. 

5.1 Results for the mathematical model 

The model presented in section [3] was implemented and used to solve the 
instances generated with 10% infeasibility rate and a time factor of 2, that 
is, 225 instances. The results were obtained using IBM ILOG Cplex 12.4 
with a time limit of 1800 seconds. 
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The model proved to be very difficult to solve. After half an hour of 
execution, it could prove optimality for only 20 instances. It also managed 
to prove that there was no feasible solution using more than one assembly 
line for 73 instances. For the remaining instances it found feasible solutions 
but with large optimality gaps. It was unable to find solutions whose cycle 
time was lower than the best known cycle time for the respective serial 
ALWABP. 

These results show that the model is inefficient in what concerns execu- 
tion time and solution quality. Therefore, we did not perform tests with the 
remaining more difficult instances. 

5.2 Results for the heuristics 

Table [2] presents the results obtained by the heuristic presented in section 
|4] for the set of instances generated, using two parallel lines. The results 
are grouped by task times variation, percentage of incompatible tasks and 
type of precedence graph. Columns C% represent the average difference 
between the best solution found by our heuristic and the best known solution 
for the serial ALWABP. Columns T indicate the average execution time of 
the heuristic in seconds. Columns P% indicate the percentage of solutions 
using parallel lines found, that is, the percentage of instances solved by 
the method. Column Best C% represents the largest improvement over the 
best know solution for the serial ALWABP among the instances with similar 
characteristics. 



Time factor 


2 


Infeasibility rate 


10% 


20% 




C% 


T(s) 


P% 


Best C% 


C% 


T(s) 


P% 


Best C% 


Bottleneck 
Chain 
Mixed 


-0,25% 
1,93% 
0,36% 


34,81 
24,47 
35,32 


58,11% 
68,92% 
64,86% 


-15,37% 
-17,56% 
-15,10% 


-2,02% 

0,23% 

-1,08% 


15,88 
14,67 
21,04 


39,19% 
45,95% 
45,95% 


-14,94% 
-15,73% 
-16,15% 


Time factor 


5 


Infeasibihty rate 


10% 


20% 




C% 


T(s) 


P% 


Best C% 


C% 


T{s) 


P% 


Best C% 


Bottleneck 
Chain 
Mixed 


8,67% 
10,09% 
11,19% 


44,88 
31,29 
47,72 


58,11% 
68,92% 
64,86% 


-21,13% 

-15,85% 
-10,67% 


5,95% 
8,64% 
8,27% 


21,48 
19,25 
29,36 


39,19% 
45,95% 
45,95% 


-10,74% 

-13,10% 

-9,96% 



Table 2: Test results for the heuristic for the ALWABP with parallel lines 

The heuristic method ran in low computational time and found better 
solutions on average, with respect to the linear model. Most of the execution 
time corresponds to the second phase of the method, that is, the evaluation 
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of the sets of workers. Moreover, the method was able to identify infeasible 
instances in less than a second. 

As it was expected, the more incompatibilities each worker has, the 
harder it is to find solutions with parallel lines. However, in many in- 
stances with low infeasibility rates, the method found solutions with lower 
cycle times than the respective solution for the serial ALWABP. Indeed, the 
heuristic found solutions with cycle time up to 21,13% lower than the cycle 
time for the best known solution for the serial ALWABP. Instances with 
bottleneck precedence graphs resulted in the lowest average difference be- 
tween the solution found by the heuristic and the respective ALWABP best 
known solution. 

Table [3] presents the results considering only improved solutions, that 
is, solutions with lower cycle time than the best known solution for the 
serial ALWABP. Column P+% shows the percentage of improved solutions. 
Column C+% presents the average improvement of the improved solutions. 
These results show that a large percentage of instances can be improved and 
that this improvement can be significative with respect to the solutions of the 
best known (often optimal) serial solutions. Average gains in productivity 
for these cases vary from 2,98% to 8,33%, depending on the precedence graph 
and worker characteristics. These results indicate that such configurations 
should be considered when planning such assembly lines. 



Time factor 


2 


Infeasibility rate 


10% 


20% 




P+% 


C+% 


P+% 


C+% 


Bottleneck 
Chain 
Mixed 


29,73% 
21,62% 
24,32% 


-6,39% 
-7,11% 
-8,23% 


22,97% 
20,27% 
22,97% 


-6,58% 
-5,98% 
-8,33% 


Time factor 


5 


Infeasibility rate 


10% 


20% 




P+% 


C+% 


P+% 


C+% 


Bottleneck 
Chain 
Mixed 


14,86% 

6,76% 

10,81% 


-5,43% 
-5,84% 
-2,98% 


10,81% 

8,11% 

12,16% 


-3,65% 
-7,22% 
-3,15% 



Table 3: Test results for the heuristic for the ALWABP with parallel lines 
(improved solutions) 
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6 Conclusions 

We have defined, modeled and solved a new assembly line balancing problem 
named Parallel Assembly Line Worker Assignment and Balancing Problem. 
This problem is motivated by the context found in sheltered work centers for 
the disabled where workers with very different characteristics are assigned 
to assembly lines. The proposed model and algorithms indicate that the 
use of parallel assembly lines in this context might improve productivity 
levels. These results encourage the study of other alternative layouts or 
the proposal of even more efficient algorithms for the planning of parallel 
assembly lines in such situations. From a practical point of view, we provide 
an additional powerful approach for coping with workers heterogeneity while 
ensuring productivity; what becomes crucial for many SWDs survival in the 
current economical context. Exploring the potential use of this approach to 
minimize the adverse effects of absenteeism is foreseen as another interesting 
research line. 
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